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Abstract 

We have made a new time-dependent calculation of the supernova production ratio of the long-lived 
isomeric state ^^'^"^Ta. Such a time-dependent solution is crucial for understanding the production 
and survival of this isotope. We include the explicit linking between the isomer and all known excited 
states. We have also calculated the properties of possible links to a conjectured excited state which 
might decrease the final isomer residual ratio. We find that the explicit time evolution of the synthesis 
of isoTa using the available nuclear data avoids the overproduction relative to ^'^®La for a u process 
neutrino temperature of 4 MeV. 

PACS numbers: 26.30.-k; 26.30.Jk; 25.30.Pt 
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The nucleosynthesis of ^^"Ta has remained an unsolved problem. For the most part this 
nucleus is bypassed by the major nucleosynthesis mechanisms of the s and r processes. This 
accounts for why this isotope is the rarest in Nature. For this reason, a variety of more exotic 
processes, such as a weak branch through excited states in ^^''Hf [1], the /3 decay of ^^^Hf 
followed by neutron capture Q], and the 7 process 3-6], have been proposed and studied 
experimentally E-^- Perhaps, the most popular scenario in recent times is ^^°Ta production in 



the u process 



nj via the ^^^Ta{ 



v,v n 



180 r 



Ta and Hf(z/e,e) Ta neutrino reactions in core- 



collapse supernovae. It is currently believed that only two isotopes (^^^La and ^^°Ta) among 



the heavy elements may be predominantly synthesized by the v process ll|]. Although the 
new calculated result based upon results of a recent nuclear experiment [l2| can reproduce the 
Solar abundance of ^^^La with charged current reactions and an electron neutrino temperature 
of 4 MeV, it overproduces the abundance of ^^"Ta. Here we investigate the possibility that 
this overestimate originates from the unique feature that the naturally occurring abundance of 
^^'^Ta is actually a meta-stable isomer (half-life of > 10^^ yr), while the true ground state is 
a 1"*" unstable state which /3-decays with a half-life of only 8.15 hr (see Fig. [1]). Therefore, a 
crucial ingredient for all of the possible production scenarios is the ratio of the population of 
the meta-stable isomer to the total production of this isotope. 

In the V process, low-spin excited states in ^^°Ta are strongly populated from ^^°Hf by 
Gamow- Teller transitions and subsequently decay preferentially to the l"*" ground state [11] (a 
similar situation for ^^^La is discussed in Ref. 13|). However, in a high temperature photon 
bath, the meta-stable isomer is excited from the ground state by (7,7') reactions through highly 
excited states. Moreover, the transition rate between the ground state and the isomer is affected 
by the changing temperature. Therefore, the final isomeric branching ratio should be evaluated 
by a time-dependent calculation. 

Previous studies have noted 

HQ 

that the observed ^^°™Ta abundance can not be in- 
ferred from their calculations until the branching between the long-lived isomer and the ground 
state is known. If one knew sufficient information on transitions between the two states of 
^^°Ta, the branching ratio could be calculated by the mathematical method of [l^. However, 
although many experiments have been carried out to identify the paths between the two states 



m 



180 



Ta 



15 



19|, such paths have never been identified and only transition probaiblities have 
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been measured [8|, |20|. In this letter, therefore, we develop a new calculation method which can 
be applied even with limited information. This method is then applied to provide a realistic 
estimate of the branching in ^'^''Ta. We show that this model leads to a good agreement between 
the final calculated isomeric abundance and the observed Solar-system abundance for a broad 
range of astrophysical parameters. 

We adopt an exponential decrease T = Toexp (— t/r) as a reasonable approximation to the 
adiabatic expansion of shock heated material in supernovae jsl. Nuclei are completely thermal- 
ized in the high temperature regime so that the population ratio of any two states is simply given 
by a quotient of their Boltzman factors, i.e. nii/mj = (2Jj + l)/(2Jj + 1) exp [— (-Ej — Ej)/kT], 
where rrii denotes the population of the state of i, with spin Jj, and excitation energy, Ei. For 
Tg = 0.1—1.0 (where Tg is the temperature in units of 10^ K) all excited states lower than a 
few hundred keV are populated. After the freeze-out each excited state decays to either the 
ground state or the isomer and must be considered 22 1. 



In the transitional region, strongly connected states are only partly thermalized. For our 
purposes we can model the excited-state structure of a deformed nucleus as consisting of two 
sets of nuclear states: 1) the ground state structure, which consists of the ground state plus 
the excited states with strong transitions to the ground state; and 2) the analogous isomeric 
structure. The transition probability between states of the two structures depends upon their 
quantum number K which is the projected component of the total angular momentum along 
the nuclear symmetry axis [23]. Transitions with AK > AI are forbidden, where AI is the 
transition multipolarity. In the case of ^^°Ta the ground state and the isomer have K = 1 and 
9, respectively. Thus, the two structures can only communicate by weak linking transitions 
and each structure can be considered independently thermalized (see Fig, |2]). We can therefore 
treat these structures as two independent nuclear species. 

To construct a model for the time evolution we consider two simple cases of linking transitions 
as shown in Fig. O In the first case there is a single linking transition between states 2 and 
3. In this case, the time- dependent evolution of the population probability of the ground-state 
structure, A^'o = X]"^f/(X]"^f given by 

dNo/dt=-P^pB23No + PrAs2{l-No) , (1) 
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where A32 and B23 are the Einstein coefficients between the indicated states, p is the photon 
density of a thermal Plank distribution, and P^^"^'' is the normalized population ratio of the 
excited state, P^^"^^ = mf^"^^/^m^^™\ In the case of kT << {E2-E1), the Einstein coefficients 
are related by pBi2 = (5'2/fl'i)^2iexp [—(£'2 — Ei)/kT], where gi = (2Jj + l) is the spin statistical 
factor, and we obtain, 

dNo/dt = -{g-s/go) exp [-(E3 - Eo)/kT]P^A32No+{93/gi) exp [-{Es - E,)/kT]P^As2{l-No) . 

(2) 

A similar expression exists for the case of a single linking transition between excited states 
4 and 5 on Fig. [2l It is straightforward to extend this to the general case of multi-linking 
transitions. The time-dependence of the population probability of the ground state structure 
is then given by 

= -E^o - exp [-{E, - Eo)/kT]A,pNo + Tp^- exp [-{E, - E,)/kT]Aj,{l - No), (3) 

where, and 1 denote the ground state and the isomer, respectively, i (j) denotes levels of 
the ground state structure (or the isomer structure) and p (g) denotes the levels of any other 
structure. 

In general, excited states of deformed nuclei are characterized by collective rotational motion. 
Each excited state is a member of a rotational band, in which the electric transition probabilities 
are enhanced by 1—2 orders of magnitude. However the interband transition probabilities 
between two such rotational bands are much hindered relative to the intraband transitions. 
Therefore, in the case of deformed nuclei, the Fj corresponding to the interband transition 
rates are much smaller than Fq and one can approximate gi/giTiTo/T f» Qi/gi^i- Finally, 
inserting A = T/h, into equation ([3]), we obtain 

^ = -E^o -^^P [-(^^ - Eo)/kT]^^No + E^rexp [-(E, - E,)/kT]^^{l - No), 
at '—^ go gi n '—^ gi n 

(4) 

We have applied this general formula [Eq. to the case of ^®°Ta. The excitation energies 



and spins have been taken from an evaluated data set of Ref. 



24| . We have calculated P^ 



,g{m) 



by taking into account all known excited states up to 600 keV excitation energy. The isomer 
population ratio = m^/^m™ is 0.85—0.94 during the transitional temperature region of Tg 
= 0.44—0.62 (see below). This result indicates that the linking transitions connected directly 
with the isomer are crucial for the total transition rate between the ground state and isomeric 
structures. Note that the transition between states 2 and 3 in Fig. [2] are negligibly small since 
the spins of the excited states above the isomer are generally larger than the spin of the isomer 
(J = 9). 

Belie et al. reported on the partial widths gi/ giV.Xo/T [meV] |20| for 9 transitions connected 
with the isomer. However, transitions between the ground state and the excited states connected 
to the isomer have not been identified (see Fig. [T]). Figure [3] shows the calculated isomer 
population ratio taking into account the 9 transitions in ^^''Ta. For initial conditions we begin 
with static thermal equilibrium among states at Tg = 1.0, and we take r = 1 s for the supernova 
temperature time constant. The population of the isomer structure decreases with decreasing 
temperature as expected. In the high temperature region (Tg > 0.62) the present calculated 
result is identical with that obtained by assuming static thermal equilibrium. However, in the 
low temperature region (Tg < 0.62) the present time-dependent calculated result is significantly 
different from the thermal equilibrium result. In the freeze-out region at low temperature 
(Tg < 0.44) the two structures are completely disconnected and the isomer population ratio 



remains fixed at Pm/ (Pm+Pgs) = 0.39 ± 0.0 



The uncertainty which is evaluated from the 



experimental errors of the energy width Fj [20| is small since both the transition rates of m ^ 
gs and gs — )■ m are proportional to Fj. 

For completeness we should also remark that it has been suggested 22|] that an unobserved 
linking transition to a state at 592 keV might also exist. To estimate the transition probability, 
we plot in Fig. H] the 9 known linking transition widths to the isomer as a function of their 
excitation energy. The transition widths increase with excitation energy and can be fit with 
a simple exponential growth. This trend can be understood in terms of a mixing effect 
caused by the Coriolis interaction, but a detailed explanation of why a simple exponential 
growth reproduces the observed transition probability so well is beyond the scope of this Letter. 
Nevertheless, from a least-square fit to this trend (dashed line on Fig. H]) we can estimate that 
an unknown state at 592 (800) keV would have a width of Qi/giTi = 0.003 (0.008) [meV]. 
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The isomeric residual population ratio is then Pm/ (Pg+Pm) = 0.18 and 0.30 for hypothetical 
additional linking transitions to states at 592 keV and 800 keV, respectively. The isomer ratio is 
sensitive to the lowest energy of the linkning transitions. Note, however, that available experi- 



mental results for this nucleus have found no evidence for direct transitions 
or 0.739 MeV {2^, and that previous 7-ray spectroscopy experiments 15 



below 1.09 MeV 



20| 



191] have measured 



no semi-direct transitions lower than 1.09 MeV; this fact indicates that the lowest energy of 
the linking transitions is the known energy of 1.09 MeV. We choose, therefore, Pm/{Pg+Pm) = 
0.39 based only on the observed 9 linking transitions for the following discussion. 

In the u and 7 processes, ^^°Ta is initially synthesized in environments with Tg > 0.62 so 
that ^^°Ta is completely thermalized and the production ratio in any previous nuclear reaction 
process does not affect the final relative isomeric population. The isomer population at freeze- 
out is insensitive to the temperature time constant because of the slow rate of change of the 
population near the freeze-out. The calculated values are almost identical for temperature 
timescales in the range of r = 0.3—3 s. As we shall now show, the only remaining sensitive 
astrophysical parameter is the neutrino temperature for the u process. 

Byelikov et al. 12| presented nucleosynthesis yields and it was concluded that ^^^La was 
reproduced by neutrino charged current reactions and z/g temperatures around 4 MeV. However, 
^^°Ta was overproduced without correcting of the isomer residual ratio. In Table 1, we present 
the modified yields taking into account the isomer residual ratio. The quoted production ratios 
of ^^^La and ^^"^Ta are normalized to ^^O which is known to be predominantly produced in 
supernovae and its solar abundance is very large compared with that of the heavy elements. 
A successful nucleosynthesis mechanism must produce abundances relative to ^^O near unity. 
Otherwise, these isotopes are under (or over) produced. Both ^^^La and ^^°Ta are now produced 
in relative amounts near unity by the charged current neutrino reactions and a temperature 
of 4 MeV. Although this result is only from the 15 Mq model the core properties do not 
depend much upon progenitor mass and the initial mass function is weighted toward low-mass 
progenitors. Hence the 15 M© model should be a good representation of a proper average over 
supernova progenitors. 

We stress that, in the present calculation, the isomer ratio is almost independent of the 
astrophysical parameters such as the peak temperature, the temperature time constant, the 
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TABLE I: Nucleosynthesis production factor relative to Solar (normalized to ^^O) for various 7 or 
process conditions (column 1) based upon the 15 Mq supernova model of Ref. [ill]. In the notation 
of Ref. {1^ . + n.c. means the 7 process plus neutral-current interactions, and ++c.c. denotes that 
charged current interactions are also included at the indicated neutrino temperatures. Column 2 gives 
the total ^^^La yields and column 4 gives the total (g.s. + isomer) ^^'^Ta yields [ill . The yields 
of ^^'''"Ta (column 3) deduced in the present work are based upon our calculated isomer residual 
population fraction of 0.39. Note that for = 4 MeV ^^'^La and ^^'^'"Ta are produced with about the 
same over-production factor. 



Model 138l^ 


180m -pg^ 


^^°Ta(g.s.+isomer) 


7-proc. only 0.190 


0.234 


0.599 


+ n.c. (6 MeV) 0.280 


0.399 


1.024 


++ c.c. (4 MeV) 1.101 


1.218 


3.123 


++ c.c. (6 MeV) 2.797 


2.109 


5.409 


++ c.c. (8 MeV) 3.222 


2.881 


7.387 



supernova neutrino energy spectrum, and the explosion energy. The final result that the over- 
production problem of ^^"Ta in the u process is largely reduced is probably robust, though only 
a full set of stellar nucleosynthesis calculations that include the time-dependent depopulation of 
the isomer will answer this question quantitatively. As shown in Table HI the present result con- 
strains the neutrino energy spectrum, which is of importance for understanding the supernova 
explosion mechanim and the detection of supernova neutrinos in near future. A previous study 



261] has shown that light element yields of the u process depend upon the neutrino oscillation 
parameter B13. Heger et al. suggested that the yields of ^'^^La and ^®°Ta is sensitive to G13 



Using the present isomer ratio, one can now study systematically the neutrino oscillation 
effect for both of light and heavy elements. 
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V process 

Y process 




FIG. 1: Partial nuclear level scheme of Ta. The ground state 



decays via /? decay with a half-life of 



20( 1 excitation energies of states with 



8.15 h, while the isomer is a meta-stable state. The measured 
weak transitions to the isomer are indicated. These excited states should also connect to the ground 
state but their decay paths are unknown. 
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FIG. 2: Schematic illustration of the nuclear structure relevant during the transitional temperature 
region. The ground state structure (i.e. the ground state and excited states above the ground state) 
is in thermal equilibrium. The isomeric structure is also in thermal equilibrium. The ground state 
and isomeric structures are connected via the indicated linking transitions. 
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FIG. 3: Calculated isomer population ratio. The solid line denotes the time-dependent calculated 
result. The dashed line denotes the isomer ratio under the condition of thermal equilibrium. In the 
high temperature regime (Tg > 0.62) both are identical. In the low temperature regime (Tg < 0.62) 
the calculated ratio is much higher than the thermal equilibrium ratio. 
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FIG. 4: Observed transition widths Qi/giTiTQ/T [meV] of the 9 known linking transitions as a function 



of their excitation energy. The Data are taken from Ref. 



20( 1 . The dashed hne is a simple exponential 



least square fit to these widths used to estimate widths for other hypothetical linking transitions. 
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